Direct measurement of topological invariants in optical lattices 
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We propose an experimental technique for classifying the topology of band structures realized in 
optical lattices, based on a generalization of topological charge pumping in quantum Hall systems 
to cold atoms in optical lattices. Time-of- flight measurement along one spatial direction combined 
with in situ detection along the transverse direction provide a direct measure of the system's Chern 
number, as we illustrate by calculations for the Hofstadter lattice. Based on an analogy with 
Wannier functions techniques of topological band theory, the method is very general and also allows 
the measurement of other topological invariants, such as the Z2 topological invariant of time-reversal 
symmetric insulators. 
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Recent advances in experimental techniques have led 
to realization of synthetic gauge fields [THS] and spin- 
orbit coupling [6l |7] in cold atomic gases. These new 
developments allow one to study a variety of topologi- 
cal phases of condensed matter physics using cold gases 
of neutral atoms trapped in optical lattices. Such topo- 
logical phases occur in systems whose Hilbert space has 
a non-trivial topological structure, and they are classi- 
fied according to the value of a corresponding quantum 
number, the topological invariant. 

Particular examples of such phases in condensed mat- 
ter include the quantum Hall insulators |8J or the quan- 
tum anomalous Hall insulators [9], where the topological 
invariant of the Hilbert space is the so-called Chern num- 
ber [10]. Some recent experiments [3-5 , llj point towards 
the possibility of realizing an optical lattice with non-zero 
Chern number in the near future. 

Once the desired lattice is created the question of 
experimental verification of the non-trivial topology 
arises. Unlike condensed matter systems, where a rou- 
tine measurement of the Hall conductance reveals the 
Chern number value [8j, cold atom systems require a 
special setup [12] for transport measurements. Edge 
states probes of condensed matter experiments also be- 
come cumbersome in a cold atom environment since the 
smooth harmonic potential washes out the edge states 
associated with the quantum Hall state. This problem 
can potentially be circumvented by enhancing a weak 
Bragg signal from the topological edge states by means 
of specifically tuned Raman transitions [T3l [T4] , Other 
approaches, which do not have an immediate analogy 
among solid state experiments, might also allow for the 
measurement of Chern numbers in optical lattices [T5l - 
[T9] . However, the quest for a universal method to obtain 
Chern numbers and other topological invariants directly 
in a single measurement, avoiding a sophisticated exper- 
imental setup or data analysis, remains open. 

In this Letter we tackle the problem of detecting topo- 
logical invariants in optical lattice systems from a very 



different perspective, drawing an analogy to the theory 
of electric polarization of crystalline solids [20j to suggest 
a simple and effective method to measure Chern num- 
bers in cold atom systems. We introduce the concept of 
hybrid time-of- flight (HTOF) images, referring to an in 
situ measurement of the cloud's density in one direction 
during free expansion in the other. The HTOF reveals 
the topology of the optical lattice just as hybrid Wannier 
functions (HWF) do in band theory [2TJ [22]. We illus- 
trate our approach by numerical simulations of a square 
optical lattice that realizes a Hofstadter model [23J, and 
discuss how it works in lattices with a more complicated 
geometry. Our method does not require the presence of 
the sharp edge states and is not affected by a soft har- 
monic trap. It can also be used to detect the Z2 topolog- 
ical invariant of time-reversal (T) symmetric topological 
insulators. 

The modern theory of electric polarization of crys- 
talline solids [20j [24] relates the electronic polariza- 
tion P to the geometry of the underlying band struc- 
ture. For a ID insulator with a single occupied band 
P=^ § BZ A(k)dk, where k is the crystal momentum, 
A(k)=i(uk\dk\uk) is the Berry connection [25] and the 
us are the lattice-periodic parts of the Bloch functions. 
Alternatively, we will use the fact that the polarization 
can be written [20J as the center of mass of the Wannier 
function constructed for the occupied band [26J, which 
can be defined as an expectation value of position oper- 
ator projected onto the occupied state [27} [28 ] . 

In two dimensions (2D) an insulating Hamiltonian can 
be viewed as a fictitious ID insulator subject to an exter- 
nal parameter k x . Polarization of this ID insulator can 
be defined by means of HWF [29J , which are localized in 
only one direction retaining Bloch character in the other. 
The polarization at each k x is given by the center of mass 
of the corresponding HWF [22] . 

The definitions of electronic polarization given above 
are gauge dependent, meaning that P is defined only 
modulo a lattice vector. For measurements one has to 
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FIG. 1. Square lattice illustrating the Hofstadter model of 
Eq. Q. The oval marks q—1 sites of the unit cell. Small ar- 
rows indicate the directions in which the phase of the hopping 
amplitude is chosen to be positive. Different colors of these 
arrows correspond to different values of the phase. Large ar- 
rows indicate the direction of the free expansion of the atomic 
cloud in hybrid time of flight (HTOF) images. 



consider the change in polarization induced by a change 
in an external parameter [20J . In the 2D insulator consid- 
ered above, k x plays the role of such a parameter. When 
k x is adiabatically changed by 27r, the change in polar- 
ization, i.e. the shift of the HWF center, is proportional 
to the Chern number [20J. This is a manifestation of 
topological charge pumping [30j I5T] . with k x being the 
adiabatic pumping parameter. 

A generalization of these ideas to cold atomic gases is 
a natural way to measure the Chern number in optical 
lattices. We replace the HWFs of band theory with hy- 
brid densities p(k x ,y), which are the particle densities 
resolved along the ^-direction as a function of k x . Note 
that, while for an extended system the HWF charge cen- 
ter position can not be reconstructed from the single par- 
ticle density [21], this becomes possible in a finite system 
where the position operator is well defined [32] . Hence in 
a system with finite extent L y in the ^-direction we can 
calculate the HWF center as 



y(k x ) 



Io V yp(k x ,y)dy 
Io V p(k x ,y)dy 



(1) 



The proportionality between the shift of the HWF cen- 
ter and the Chern number still holds in the open system. 
The Chern number measures the charge transported from 
one boundary to the other as k x is cycled by 2tt. An ex- 
perimental measurement of the shift in the hybrid density 
will hence directly determine the Chern number. 

Experimentally, p(k Xl y) is measured by an HTOF 
measurement, in which the lattice and trap are switched 
off along the ^-direction, while keeping the lattice and 
harmonic confinement unchanged in the ^-direction. In 
the long time limit [33J TOF images map out the crys- 
tal momentum distribution along k x [34J. At the same 
time, the system is still confined in the ^-direction and a 
real-space density resolution can be done. 

We now show HTOF unambiguously determine the 
Chern number by performing a numerical simulation of 




FIG. 2. Results for a ribbon (L y — 10) of the Hofstadter lat- 
tice model with p/q = 1/7 and corresponding Chern number 
C = 1. (a) Energy spectrum: two edge states (in red) cross 
the bulk energy gap. (b) Hybrid density p(k x , n) shifts by one 
unit cell in a 27r-cycle. (c) The center of mass of the hybrid 
density Eq. (|3| jumps by one unit cell. 



the Hofstadter model [23J on a square lattice. Its Hamil- 
tonian is given by 



lattice 



= -J, 



m,n+l^rn,n 



+H.c, 



(2) 

where J a is the hopping amplitude in the a={x, y} di- 
rection and c m ^ n is the fermionic annihilation operator, 
with m and n being the column and row indices of the 
lattice (see Fig [I]). 

An atom hopping clock wise around a plaquette accu- 
mulates a phase <£. We consider & = p/q where p and q 
are two relatively prime integers. The hopping matrix el- 
ements J x e l27rn ® along the x-direction depend on the row 
index n so that each unit cell contains q sites. In the fol- 
lowing we fix q=7 and assume that only the lowest band 
is occupied. The Chern number C of the lowest band 
is determined by the Diophantine equation l=qs+pC 
[10, 35, 36J, where s is an integer and \C\ < q/2. 

We first consider an infinite ribbon of this model with 
width L y =10 setting J x =J y =J, and p = 1 which cor- 
responds to a Chern number C = 1. In the spectrum 
shown in Fig. [2ja) we see that, as expected for p = 1, 
two edge states cross the Fermi level. Analogously to the 
2D insulator considered above, this setup can be viewed 
as a finite ID chain subject to a /^-driven pump. From 
this point of view the hybrid density p(k x ,n) describes 
the change in the density of the ID system as a function 
of the pumping parameter. Figure [2^b) shows that the 
hybrid density is shifted by one unit cell in the bulk, indi- 
cating that a single charge is pumped across the system, 
as expected for C = 1. 

We calculate the k x dependence of the HWF center by 
taking a tight-binding limit of Eq. Q 

2l n p{kx,n) 
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As shown in Fig. [2jc), n(k x ) jumps by one unit cell 
(q=7 sites), analogous to the HWF shift in Chern in- 
sulators [21J. 

Having established a clear connection between the hy- 
brid density and Chern number, we now turn to a more 
realistic case by adding a harmonic trapping potential of 
the form 



Hi 



trap 



V T V[(m - L x /2) 2 + (n - qL y /2) 2 ]J 



m,n 



where L q 



(y) 



(4) 

denotes the number of unit cells in the 
x(y) direction. The system contains qL y rows and L x 
columns. The values of Vr and the number of atoms 
N — 300 are chosen such that the atom cloud has a large 
insulating region corresponding to 1/q filling at the trap 
center. We now consider the cases p=l and p=3, corre- 
sponding to Chern numbers +1 and —2 respectively. 

There is no significant difference in the real space den- 
sities of the two states since the harmonic trap smears 
out the edge states |13| . In contrast, the HTOF den- 
sity profiles allow one to directly read of the Chern num- 
bers. We calculate these HTOF images by first solving 
the Schrodinger equation (lattice + H trap )^i(m,n) = 
£ii/ji(m,n) and constructing HWF's by means of the 
Fourier transform in the ^-direction 

L x -1 

*l>i(k x ,n) = J2 e lkxm Mm,n). (5) 

ra=0 

We then use these wavefunctions to construct the hybrid 
particle density of the HTOF measurement 



1 N 



(6) 



HTOF images obtained in this way are shown in Fig. [3] 
and clearly exhibit the topological charge pumping effect 
despite the presence of a trap: the hybrid density shifts 
by C unit cells along the ^-direction as k x changes from 
— 7r to 7r. Thus, the hybrid density is an accurate probe 
of topological properties and allows to directly measure 
the Chern number. 

To get deeper understanding we consider the case of 
vanishing transverse coupling 7^=0, corresponding to a 
set of decoupled tight-binding chains with dispersions 
£n(kx) = — 2 J x cos(k x — 27rnp/q). The position of the 
band minimum shifts by 2irp/q from one chain to the 
next, as shown in Fig. [4^a) for p=3. If an infinitesi- 
mal coupling J y is added, the 2D lattice is in the Chern 
insulating regime. Charge pumping can be inferred by 
tracing the change in the position of the valence band 
minima for weakly coupled chains: connecting nearest 
neighbor points in Fig. |4^a) reveals a shift by two unit 
cells (\C\=2) in the course of a 2tt change of momentum 
(dashed red arrow in Fig. [4ja)). This illustrates the geo- 
metrical interpretation of the Diophantine equation [10] , 
as also discussed in Ref. [37] in a different context. 





FIG. 3. The HTOF images for the Hoftadter lattice (L x = 70, 
L y = 10, N = 300) in the presence of the harmonic trap: (a) 
for p/q = 1/7 with C= 1 (V T = 0.001), (b) for p/q = 3/7 with 
C = — 2 (Vt = 0.001). Chern numbers can be determined as 
the number of unit cells traversed by the hybrid density in 
the course of a 27r-cycle. Upward (downward) direction of the 
shift corresponds to positive (negative) Chern number. The 
broadening of lines correspond to exponential localization of 
the peaks of the hybrid density. 



The above analysis allows for an alternative way of 
describing the HTOF results presented in Fig. |3ja-b). 
We introduce sublattice densities p a {k Xl h) 1 which corre- 
spond to the particle density on the a-th site of the n-th 
unit cell. These sublattice densities, shown in Fig. |4^b), 
shift along the ^-direction as k x changes, illustrating the 
motion of charge. The motion of charge can also be 
tracked in the total sublattice density obtained by sum- 
ming p a {k Xl h) over the unit cells: 



Ly-l 

n=0 



(7) 



Thus, the topological nature of the state can also be seen 
in the total sublattice density, which is potentially acces- 
sible in a TOF experiment that can distinguish different 
sublattices |3, 38J. However, such an analysis is specific 
to the Hofstadter model, while the HTOF measurement 
is generally applicable. 

The square lattice considered so far is particularly sim- 
ple, since a straightforward partition into ID chains is 
possible. HTOF measurement can also work for other 
lattice geometries such as the honeycomb lattice, which 
is topologically equivalent to the brick- wall lattice shown 
in Fig. (5ja). Such a lattice was used to create Dirac 
points in optical lattices [11]. This lattice is a potential 
candidate for realizing the Haldane model [9j, which is a 
canonical example of a Chern insulator. 

A partition of the brick-wall lattice into ID chains, 
along which the charge is pumped, is illustrated in 
Fig. [5|a) with solid dark bonds. The chains consist of 
two sublattices offset from each other in the x-direction. 
Due to this offset the charge pumping is not directly visi- 
ble in HTOF image unless one separately measures them 
for each sublattice (for example using the superlattice 
technique of Ref. |38j). As illustrated in Fig. the 
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FIG. 4. (a) Valence band minima of decoupled (J y = 0) ID 
chains of the Hofstadter model with p/q = 3/7. Different 
colors correspond to different sublattice chains. Solid grey 
arrows indicate the shift in the dispersion due to the phase 
factors. Dashed red arrows connect nearest neighbor points 
and illustrate charge pumping, (b) The sublattice densities 
p a (k x ,n) within a unit cell and the total sublattice density 
M a {k x ) for different values of a. The color scheme is the 
same as in left panel and the model parameters are the same 
in Fig.^b). 
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FIG. 5. (a) The brick- wall lattice split into ID chains. Thick 
(black) and light (yellow) bonds indicate intra- and inter- 
chain hoppings respectively. The shape of an imposed super- 
lattice potential Vs(x) is shown schematically (in gray), (b) 
Energy spectrum of a Haldane model in the ribbon geometry. 
Edge states are shown in red. (c) The shift of the center of 
mass of the hybrid density indicates non-trivial Chern number 

(c = i). 



center of mass of the hybrid density along the zig-zag 
ID chain is indeed shifted by one unit cell along the ID 
chain, revealing the Chern number C = 1 of the Haldane 
model. 

The HTOF technique can not only determine Chern 
numbers, but also the Z2 topological invariant of T- 
symmetric insulators [39J. In band theory this invariant 
can be obtained by means of HWFs [40-42J, once they 



form T images of one another. The occupied single par- 
ticle states of the insulator can always be split into two 
sets of states (related by T symmetry) where each set has 
a well defined Chern number [43J. For a wide range of 
models such a splitting can be achieved by projecting the 
occupied states onto particular spin directions [44 , 45J. If 
the values of thus obtained spin Chern numbers are odd, 
the system is in the Z2-insulating regime. In the context 
of cold atoms, considering a minimal model with only 
two occupied bands, a spin-projected HTOF measure- 
ment of spin-resolved densities would serve as a direct 
measurement of the spin Chern numbers, and hence of 
the Z2 -invariant. 

The proposed HTOF technique is not only practical, 
providing exhaustive information about the topological 
state of an optical lattice, but is also a conceptually novel 
idea for using a cold atom lattice as a quantum simula- 
tor. Hybrid density measurement as proposed here for 
cold atom systems are not possible in condensed matter 
experiments and the HWFs are used only in computer 
simulations. These experiments can thus access novel 
probes of topological order and will give rise to further 
implementations of so far numerical experiments of con- 
densed matter in real experiments on cold atom systems. 
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